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Abstract. We review the possible roles of large scale shocks as particle accelerators in 
clusters of galaxies. Recent observational and theoretical work has suggested that high 
energy charged particles may constitute a substantial pressure component in clusters. 
If true that would alter the expected dynamical evolution of clusters and increase the 
dynamical masses consistent with hydrostatic equilibrium. Moderately strong shocks 
are probably common in clusters, through the actions of several agents. The most 
obvious of these agents include winds from galaxies undergoing intense episodes of 
starbursts, active galaxies and cosmic inflows, such as accretion and cluster mergers. 
We describe our own work derived from simulations of large scale structure formation, 
in which we have, for the first time, explicitly included passive components of high 
energy particles. We find, indeed that shocks associated with these large scale flows can 
lead to nonthermal particle pressures big enough to influence cluster dynamics. These 
same simulations allow us also to compute nonthermal emissions from the clusters. 
Here we present resulting predictions of 7-ray fluxes. 



INTRODUCTION 



Clusters of galaxies are very important probes for cosmology, since they are the 
largest bound systems in the universe. They represent the nonlinear development 
of large-scale perturbations in the early Hubble flow. Matter, energy and entropy 
can be added from outside through infall or cluster mergers, but with the excep- 
tion of a small fraction of the energy as photons, neutrinos and ultra high energy 
cosmic-rays, nothing leaves clusters once they form. Thus, they provide unique 
records of the history of the universe. The statistics of cluster masses and their 
dynamical properties, including, for instance, the relative proportions of baryonic 
and nonbaryonic matter, are commonly used to test basic cosmological models [1]. 
While galaxies are the most obvious constituents of clusters in visible light, most 



of the matter in clusters is nonbaryonic, and even the baryonic matter is primarily 
contained within the diffuse intracluster medium (ICM), rather than the galaxies. 
The temperature and density distribution of the ICM gas directly reflects the dy- 
namical state of the clusters, as well as holding a history of important dissipative 
processes encountered by the gas [39]. Elemental abundances in the ICM are seen 
as key indicators of the star formation histories and of galaxy evolution, more gen- 
erally [43,53]. The dynamical states of clusters have also received much attention. 
While cluster ICMs sometimes appear relaxed, often the situation is quite differ- 
ent, with clear indications of high speed flows [13,15,32,33,40] demonstrating that 
cluster environments can be violent. 

The possible importance of nonthermal components in the ICM has recently 
raised great interest. There is growing evidence that magnetic fields [9,35], and 
high energy charged particles [6,27,38] may constitute significant dynamical com- 
ponents of the ICM in at least some clusters. If these components are generally 
strong, they would impact on a wide range of issues, beginning with estimates of 
cluster masses derived from assumptions of ICM hydrostatic equilibrium. In addi- 
tion, since high energy charged particles and magnetic fields do not readily radiate 
away their pressures, they could tend to inhibit cooling flows. Radio, EUV and hard 
X-ray emissions resulting from energetic electrons have already stimulated much 
discussion about cluster physics as well as the evolution of the clusters and their 
galaxies [4,19,34,45,49,53]. While high energy 7-rays have not yet been detected 
from clusters, recent estimates of 7-ray luminosities from high energy particle in- 
teractions in the nearest rich clusters, such as Coma, are within the range of what 
may be detected in the next generation of 7-ray observatories [6,10]. 

Our focus here is the possible relationships between shock waves generated by 
high speed flows in clusters and the origins of the high energy charged particles. 
For convenience, we call all such particles cosmic-rays or "CRs". A few basic 
considerations of the CRs are necessary to set the stage. So far direct observations 
reveal only the presence of CR electrons, although one generally assumes that 
protons are present with at least comparable numbers and greater energy content 
[38]. Relativistic protons below energies where photo-pion production through 
interaction with the cosmic microwave background (CMB) becomes important, 
(~ 10 9 - 5 GeV) do not lose significant fractions of their energy in a cluster over 
a Hubble time [4]. Also, up to somewhat lower energies they should scatter off 
magnetic irregularities sufficiently to remain trapped in the cluster for a similar 
duration [53]. Relativistic electrons, on the other hand can only survive Coulomb 
losses and radiative losses for more than ~ 1 Gyr in a relatively narrow energy 
band centered roughly around 50 MeV [50]. Thus, observed relativistic electrons 
must be freshly accelerated, reaccelerated or be introduced as decay products of 
other interactions [8,12,16]. 

The Coma cluster provides a convenient scale for the energetics involved in pro- 
ducing and containing the CRs. The dynamical mass estimates of Coma are a bit 
over 10 15 Mq [18], leading to an ICM mass ~ 10 14 Mq. With a mean ICM temper- 
ature of 8 keV [55] the total thermal energy in the ICM ~ 4 x 10 63 erg. Observed 



excess EUV emission from Coma, interpreted as inverse Compton scattered CMB 
photons, leads to an estimated ^ 10 61 erg in total CR electron energy and, through 
conventional arguments, > 10 63 erg in protons [38,49]. If we accept that figure for 
argument and assume the CR protons were supplied continuously over a Hubble 
time the mean CR energy input rate would be ~ 3 x 10 45 erg/sec. By comparison 
the energy input to CRs in our own galaxy is only about 10 41 erg/sec [14,17]. This 
estimate of CR energy in Coma makes it clearly very important dynamically, and 
the implied rate of energy input severely restricts the range of possible accelerators. 



LARGE-SCALE SHOCKS AS ACCELERATORS 

There are likely to be many particle accelerators in clusters, but only a few seem 
really capable of accounting for the CR energy content implied above for clusters 
like Coma. Thus, for example, while second order Fermi acceleration from MHD 
turbulence in the ICM might conceivably play a role in electron acceleration [8,16], 
that would not seem to be a very fruitful approach to accounting for proton CRs, 
if their energy content exceeds ~ 10 63 erg in rich clusters. In fact, most current 
models for proton acceleration depend on first order Fermi acceleration resulting 
from CR diffusion near shocks. This so-called "diffusive shock acceleration", or 
DSA, can be efficient enough to put into CRs several tens of percent of the total 
energy flux through a shock [3] . Following that lead, the initial task becomes one of 
identifying shocks that dissipate enough energy to explain the energetics suggested 
for Coma and by implication other rich clusters. It turns out that large scale 
shocks are probably fairly common in ICMs, so it may indeed be reasonable to 
expect a substantial energy in CR protons there. Electrons are more difficult to 
accelerate and maintain than protons. First, they radiate energy much more rapidly 
[50]. In addition, they are harder to get started in most accelerators, because of 
their small gyroradii at suprathermal energies. Still, from galactic CRs and other 
environments, such as radio galaxies, we have abundant evidence that electrons 
can be accelerated. Also, inelastic proton-proton interactions with the ICM can 
generate a significant secondary electron-positron component through pion decay 
[11,12]- 

Several candidate classes of large shock structures come to mind immediately 
in clusters. For example, the facts that ICM gas is generally enriched with nucle- 
osynthetic products to 10 — 30% of the solar value [43] and that starburst galaxies 
like M82 produce strong winds, have been cited to support the idea that termina- 
tion shocks in galactic winds may have been common enough during early "bright 
phases" of galaxy evolution to accelerate a substantial high energy CR population 
[26,53]. Starburst-like energy deposition rates are probably necessary to generate 
winds inside rich clusters, in order to overcome the substantial thermal pressure of 
the ICM [53]. If one takes an ICM Fe mass ~ few x 10 11 Mq [7] and assumes it 
resulted from supernovae, then ~ 10 12 events, corresponding to about 10 63 erg of 



kinetic energy would have been released into a rich cluster like Coma [54]. That 
is close to the CR energy content required, but does not yet account for a realistic 
conversion efficiency Even if we assume ~100% of the supernova energy ends up 
as kinetic energy in the winds and that we can carry over from galactic supernova 




FIGURE 1. Properties of a simulated Mach 5 shock including backreaction from nonthermal 
particles accelerated by DSA after being injected through "thermal leakage". The momentum 
distribution of the particles just downstream of the shock is also shown. The dotted lines show 
initial conditions, while the solid/dashed lines show properties after CR protons have been ac- 
celerated to energies ~ 100 GeV and the shock structure is almost steady. Bohm-like spatial 
diffusion, k oc p, is assumed. 



a CR acceleration efficiency ~ 10 — 20% [22], this model seems to struggle to 
account for the very large energy contents if the recent EUV observations represent 
inverse Compton emission and the implied CRs fill the entire cluster [38]. On the 
other hand, these winds may still be effective accelerators or reaccelerators at a 
somewhat lower level, as well as contributors to the ICM magnetic field [34]. 

CRs escaping from radio galaxies have long been seen as potential sources of 
relativistic particle populations in clusters and sometimes beyond [19,23]. These 
CRs are usually seen to be accelerated by the terminal shocks of the jets powering 
the radio lobes, but it also seems plausible that the bow shocks of at least the 
most powerful radio jets in cluster environments might be effective accelerators. 
Energetically, radio galaxies come close to required input levels if we can assume 
on average that a rich cluster contains at least one at all times [4] . The most lumi- 
nous radio galaxies are estimated to be driven by jet kinetic powers ^ 10 45 erg/sec 
[2], which is comparable to our earlier estimated average rate of CR production 
for Coma. Individual radio galaxy lifetimes are generally estimated at ~ 10 8 years 
[2], so such a cluster would need as many as ~ 10 2 luminous galaxies capable of 
generating a powerful radio jet or a fairly high duty cycle in a smaller number of 
galaxies. On the other hand, this localized source model may have some difficulty 
explaining the diffuse radio halos seen in some luminous clusters [37], including 
Coma, because of large energy loss rates for high energy electrons and because low 
energy electrons probably cannot diffuse very far from their source [23]. Reaccel- 
eration can, of course, reduce or eliminate this issue [8]. Similarly, if the energy 
were deposited in high energy protons the difficulty is reduced. Then the nonther- 
mal emissions observed so far would presumably be explained through secondary 
electrons. There are also suggestions that the relatively light plasma in radio lobes 
may remain confined in buoyant "ghost" bubbles, so that they can be reenergized 
and illuminated at later times by passage of other large scale shocks [20]. 

Another scenario gaining some strong recent support is CR acceleration at very 
large and long-lived shocks resulting from cosmic structure formation. The likely 
existence of strong accretion shocks several megaparsecs from cluster cores has 
been recognized for a long time [5,47], and they have been suggested as possible 
sources of ultra-high energy CRs [29,30] as well as seeds for the ICM magnetic 
field [36]. Shocks resulting from discrete cluster merger events have also been 
recognized through X-ray structure in clusters [13,40], and cited as particle accel- 
erators to account for diffuse cluster radio halos [21,52] or so-called "relic" radio 
sources [45,44]. The overall energetics of these shocks is generally acceptable for 
the production of the CRs needed. Typical flow speeds in and around clusters 
will be vj ~ (2G M c i/ Rd) 1 / 2 ~ 2 x 10 3 km/sec, leading to an available power 
~ pciv)R\ ~ 10 46 erg/sec, using M d ~ 1O 15 M and R d ~ 2 Mpc. Accretion 
shocks far from cluster cores can be of very high Mach number and are responsi- 
ble for initial heating of the ICM. They appear, however, to be less important as 
potential sources of CRs inside clusters than weaker, "internal" shocks associated 
with mergers and other flows penetrating deeper into the clusters. The reason is 
that the latter shocks repeatedly process the ICM material, whereas the accretion 



shocks do it only once. 

Miniati et al [41] explored the statistics of shock heating in both SCDM and 
ACDM cosmologies. They found that far more mass gets processed through the 
internal shocks than the accretion shocks. Their analysis showed that the most 
common shock encounters by the ICM material involve shock Mach numbers less 
than 10, with the peak around M ~ 5. That is significant, since such shocks are 
strong enough to transfer as much as 20 - 30% of the postshock pressure to CRs. On 
the other hand such shocks are only modestly modified in their structures through 
backreaction of CRs diffusing upstream. Thus, so-called "test particle" estimates 
of the CR distributions are a reasonable first guess to what we may expect. Figure 
1 illustrates the evolution of a CR modified Mach 5 shock as computed from a fully 
nonlinear DSA simulation carried out using methods detailed elsewhere [31]. 

SIMULATING CR ACCELERATION IN COSMIC 
STRUCTURE FORMATION 

So far, only a few clusters have shown direct evidence of large shocks. However, 
recent cosmic structure formation simulations have demonstrated that large scale 
shocks in clusters are probably much more common and much more complex than 
these simple perceptions suggest [41]. Since clusters tend to form at the intersec- 
tions of cosmic filaments, they accrete matter in very unsteady and nonisotropic 
patterns. In addition to discrete cluster mergers, general, larger scale flows and 
shocks associated with them propagate down the filaments and through the clus- 
ters. When cluster mergers take place the accretion shocks associated with the 
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FIGURE 2. Mean Pcn/Pth measured in simulated clusters inside r = 0.5 Mpc/h at z = 0. 



individual clusters add to the shocks that form in direct response to the merger. 
The net result of all of this is a rich web of shocks, which often penetrates into the 
inner regions of the clusters. 

The recent detections of shocks from merger events has been noted above. But 
the existence of accretion shocks is not necessarily easy to detect directly. On the 
other hand, infalling clouds of unshocked, warm gas outside clusters may be identi- 
fied through the absorption lines of quasars located inside clusters. This warm low 
density gas of 10 4 — 10 5 K is photoionized by the diffuse radiation from the hot post- 
shock gas and the diffuse cosmic background radiation. In some studies [24,56], the 
C IV absorption systems of quasar emission lines with \z abs — zq SO \ < 3000 km s" 1 
are interpreted as clouds associated with rich clusters where the quasars reside. 
The characteristics of these systems of C IV absorbers are different from those 
of the typical intervening C IV absorbers. It has been noted that the veloc- 
ity difference is unexpectedly large compared to the typical velocity dispersion 
of galaxies (400 — 1200km s _1 ) in rich clusters [24]. An accretion velocity, how- 
ever, is a bit larger than the galaxy velocity dispersions, since it is given by 
Vacc = 1-31 x 10 3 km s~ 1 (T c //6.06keV) 1 / 2 , for example, in an Einstein-de Sitter uni- 
verse [5,46]. Thus it is possible that these absorption systems are in fact the infalling 
clumps of gas upstream of the accretion shock. 

Large shocks may be very efficient sources of CRs in clusters. We now outline 
a computational study of this phenomenon as revealed in a simulation by Miniati 
[42]. Briefly, this calculation used an Eulerian "TVD" hydro + N-body cosmology 
code [48] with passive treatments of magnetic fields [36] and CRs added. The CR 
protons and electron populations were accelerated according to DSA test-particle 
theory, then evolved to include adiabatic, radiative and Coulomb energy losses up 
to 10 6 GeV, using methods similar to those described previously [25]. The upper 




FIGURE 3. Volume renderings of clusters formed in a SCDM simulation. Left: Shock surfaces. 
Right: Thermal bremsstrahlung emissivity. A 14 Mpc/h 3 volume is shown at z = from a 50 
Mpc/h 3 simulation. 



energy limit was chosen so spatial diffusion could be neglected in order to reduce 
computational costs, and since most observable emissions can be studied in this 
regime. 

Secondary electron production due to pion decays, was also included. CRs were 
injected at shocks according to the "thermal leakage" model described by Kang & 
Jones [28] . Roughly speaking, a small fraction of the downstream thermal particles, 
whose velocities are large enough for them to escape back across the shock, are 
assumed to be "injected" into the nonthermal particle population, whence they are 
subject to DSA. A 50 Mpc/h 3 periodic section of a SCDM universe was simulated 
in this case on a 256 3 grid, with h = 0.5, Vl M — 1, &b — 0.13 and <r 8 = 0.6. 
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FIGURE 4. Vertical slices through the clusters shown in Figure 1. Left: Thermal 
bremsstrahlung cmissivity with shock surfaces superposed. Right: Shock surfaces superposed 
on the projected flow velocity field. 




FIGURE 5. Vertical slices as in Figure 3. Shown are Left: Log thermal proton density. Right: 
Log CR proton density. 



Figure 2 shows the mean ratio of CR pressure to thermal pressure in the central 
regions of clusters in this simulation. These results are indicative only, since the 
CRs were not coupled back dynamically in cluster evolution. The mean ratio, 
Pcn/Pth ~ 0.25, however, consistent with the results of nonlinear Mach 5 CR shock 
simulations mentioned earlier. There is considerable scatter, since each cluster has 
a unique shock history. These nonthermal pressures are high enough that they 
would produce non-negligible influences on the internal cluster dynamics. Clearly, 
future simulations should consider those effects self-consistently. 

The shock structures associated with the clusters are very complex. Figure 3 
isolates from the simulation a 14 Mpc/h 3 cube at z = 0, centered on a pair of 
interacting clusters. It shows clearly the web of shock surfaces in and around the 
clusters. The manner in which these shocks have formed and the considerable de- 
gree to which they penetrate into the cluster cores is shown in Figure 4. Here we 
have taken a vertical slice through the clusters in Figure 3, and show the inter- 
sections of the shock surfaces (contours of V • v ) with the slice plane. These are 
superposed onto the thermal bremsstrahlung X-ray emissivity and projections of 
the vector velocity field. Associated cosmic filamentary features are evident, as well 
as the fact that strong flows are focussed along the filaments and directed at the 
clusters. Where those flows enter the clusters shocks sometimes penetrate deeply 
into the cores. 

Figure 5 provides a comparison of gas and CR distributions in the same plane. 
On the left we show the spatial distribution of thermal protons, and on the right 
the distribution of CR protons. While there is a crude correspondence between 
the CRs and the gas, they clearly are not identical. The most obvious distinction 
is the existence of sharp edges in the CR distribution, reflecting their injection at 
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FIGURE 6. High energy 7-ray flux (> 100 MeV) from simulated clusters at z = assuming a 
distance 70Mpc/h, with h = 0.5. 



shocks, while the thermal density increases smoothly towards cluster centers. The 
rate of 7-ray production from pion decay depends on the collision rate between 
thermal and CR protons, so scales as the sum of these two (logarithmic) images. 
Thus, the 7 ray emission may have extensions beyond the thermal X-rays, and may 
show fine structure not evident in the thermal emission, as well. As mentioned 
earlier, most of the shock processing of the gas in clusters involves shocks with 
Mach numbers ~ 5 or so. Thus, the computed CR energy distributions typically 
resemble N(E) oc E~ 21 , and the predicted 7-ray spectrum has a similar form. 

In Figure 6 we show the simulated 7-ray flux, -F 7 (> lOOMeV), emitted within 1.3 
Mpc/h of cluster centers as a function of the cluster temperature, T x , assuming each 
is at the distance of the Coma cluster, 70 Mpc/h with h = 0.5. The least squares 
fit to this distribution gives F 7 oc T! 2 ' 95 . A scaling close to F 1 oc makes sense 
in this model from the fact that the kinetic power in accretion shocks of virialized 
clusters scales as Tj, providing a scaling for P C r or Ucr, while the mean cluster 
baryon densities scale approximately as T x in the simulation. The F^:T X relation 
found in this calculation is considerably steeper than in some other models, such as 
that by Colafrancesco & Blasi [10], which assumes all the CR protons diffuse from 
a central injection point and predicts closer to F 7 oc T° 5 . Our simulation utilized a 
box too small to produce clusters as massive as Coma, but an extrapolation of the 
7-ray relation in Figure 6 to T x = 8.3 keV leads to a prediction of about 3 x 10~ 9 
counts/sec/cm 2 , roughly the same as the value found by Colafrancesco and Blasi 
for the same cluster mass, despite the other differences in the models. The current 
EGRET limit for Coma is 4 x 10~ 8 [51], so an order of magnitude higher. Still, 
the coming generation of 7-ray telescopes should be able to test such models as 
this. Indeed, those telescopes promise a genuine opportunity to probe the physics 
of galaxy clusters in telling ways not possible before. 
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